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Ambiguity in renormalization of X-j unction 
between Luttinger liquid wires 
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Abstract 

Wt study four-lead junction of semi-infinite wires by using fermionic 
representation in the scattering state formalism. The model of spinless 
fermions with short-range Luttinger liquid type interaction is used. We 
find the renormalization group (RG) equations for the conductances of 
the system in the first order of fermionic interaction. In contrast to 
the well-known cases of two-lead and three-lead junctions, we show that 
the RG equations cannot he expressed solely in terms of conductances. 

The appearing ambiguity is resolved by choosing a certain sign defined 
at the level of S-matrix, but hidden in conductances. The origin of 
this Z2 symmetry is traced hack to the particle-hole symmetry of our 
Hamiltonian. We demonstrate two distinct RG flows from any point 
in the space of conductances. The scaling exponents at the fixed points 
are not affected by these observations. 

Keywords: renormalization gronp, Lnttinger liqnid, junctions, electrical 
conductance 


1 Introduction 

One-dimensional (ID) systems of interacting fermions served as toy models 
for theorists since more than hfty years. |1] The advances in technology and 
fabrication during last two decades propelled these models into the forefront 
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of current research in physics at nanoscale. Examples of such systems include 
carbon nanotubes, |2] ID edge states in quantum Hall regime [3] and in topo¬ 
logical insulators. mu One dimensional conductors and junctions between 
them should be important ingredients of any future electronic devices. [6] 
It is well known, that the transparency (conductance) of such junctions to 
the electric current is subject to renormalization due to interaction between 
fermions within the wire. For physical case of electronic repulsion this con¬ 
ductance tends to zero at low temperatures and applied voltages according to 
power law, with the exponent dehned by the interaction strength. There are 
two approaches developed for the description of this conductance renormal¬ 
ization. One approach uses the so-called bosonization technique, which treats 
the interaction between electrons in the bulk of the wire exactly, and consid¬ 
ers the junction as a boundary. [7] The system is described in terms of chiral 
fermionic densities (currents), and the boundary conditions are imposed onto 
these currents. Fixed points of the boundary action are determined and per¬ 
turbations around them are classihed according to their scaling dimensions. 
m This way one judges about (in)stability of various fixed points, i.e. the 
limiting values of conductance obtained during renormalization procedure. 

Another approach to renormalization of conductance was hrst formulated 
in the limit of weak interaction in the bulk. In this approach one describes an 
arbitrary junction via its unitary S-matrix, dehned in the absence of interac¬ 
tion. in For N semi-wires met at one junction this matrix belongs to U{N) 
group. The fermion wave function are given in a scattering state formalism, 
and corrections to S-matrix are subsequently calculated. m These correc¬ 
tions are logarithmically divergent which eventually leads to the renormaliza¬ 
tion group (RG) equation for the S-matrix and conductance. This approach 
was essentially improved in [T21 [13] by taking into account higher orders of in¬ 
teraction and subleading logarithms in perturbation theory for conductance. 
By summing up the appearing series one obtains the non-pertrubative RG 
equation for the conductance, whose solutions nicely reproduce those scaling 
exponents which are known exactly from the bosonization approach. 

Up to now the main focus of theoretical analysis concerned the simpler 
cases N = 2 and N = 3, which describe an impurity in one wire and Y- 
junction with three leads, respectively. The case of four-lead junction with 
iV = 4 received less attention, partly because of the difficulties in the de¬ 
scription of A-matrix and ensuing analysis of RG equation even in the lowest 
order of perturbation theory. We mention her recent studies for a special 
model cases of R-matrix with iV = 4 while discussing regular networks of 
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Luttinger liquids [2], superconducting hybrid junctions [2], and tunneling 
between two helical edge states of the topological insulators [T6] . 

We notice that in all previous studies of two-lead and three-lead junctions 
the RG equations written in terms of S'-matrix and in terms of conductance 
matrix (dehned by squares of matrix elements IS'jjp) were equivalent. The 
same equivalence held in the considered simpler model cases of four-lead 
junctions as well. However, there is no one-to-one correspondence for > 4 
between the matrix of conductances and the S'-matrix, even after removal of 
trivial phase factors. Mathematically, it is known as an ambiguity in restoring 
the unitary matrices from the unistochastic ones. dzi ng In general, this 
ambiguity can even be continuous and three-dimensional, while in our model 
below we only hnd a double discrete ambiguity which is a proven minimum 
for symmetric S'-matrix in f/(4) group. [2] 

The main goal of our paper is to demonstrate this ambiguity already in 
the hrst order RG equations. We show for a particular physically motivated 
class of S'-matrices the description of the junction in terms of conductances 
is incomplete. We hnd two possible RG hows starting from any point in the 
conductances’ space. The RG hxed points and the scaling exponents are not 
inhuenced by this ambiguity. 

The plan of the paper is as follows. We dehne our model in Sec. we 
explain the notion of reduced conductances in Sec. the RG equations are 
discussed in Sec. The ambiguity of the RG equations is revealed and 
discussed in Sec. we present our concluding remarks in Sec. 


2 The model 

We consider a model of interacting spinless fermions describing two quantum 
wires connected by a junction in the middle of the wires. Alternatively, we 
speak of four semi-wires connected at a single spot. In the continuum limit, 
linearizing the spectrum at the Fermi energy and including forward scattering 
interaction of strength Qj in wire j, we may write the Tomonaga-Luttiger 
liquid Hamiltonian in the representation of incoming and outgoing waves in 
lead j (fermion operators '^j,out) as 
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H^nt ^ 27iVF''^gjPjPjQ{x; —L, —1). 


Here V’ 2 ,m, V’a.m,'04,m) denotes a vector operator of incom¬ 

ing fermions and the corresponding vector of ontgoing fermions is expressed 
throngh the S'-matrix as "^outix) = S ■ '^in{—x) . In the chiral representation 
we are using positions on the negative (positive) semi-axis corresponding to 
incoming (outgoing) waves. We consider quantum wires of finite length L, 
contacted by reservoirs. The transition from wire to reservoir is assumed to 
be adiabatic (i.e. produces no additional potential scattering). The junc¬ 
tion is assumed to have microscopic extension I of the order of the Fermi 
wave length. Interaction effects inside the junction are neglected. This is 
expressed by the window function ©(x; —L, —1) = 1, if —L < x < —I, and 
zero otherwise. The regions x < —L are thus regarded as reservoirs or leads 
labeled j = 1,2, 3,4. We put the Fermi velocity Ui? = 1 from now on. The 
interaction term of the Hamiltonian is expressed in terms of density opera¬ 
tors pj^in = T+pjT = pj, and pj^out = = pj, where pj = S~^ ■ pj ■ S and 

the density matrices are given by {pj)a^ = SapSaj and {pj)ai 3 = The 

S'-matrix describes the scattering at the junction and belongs to t/(4) group. 

3 Reduced conductances 

In the linear response regime our system is characterized by the matrix of 
conductances defined by Jj = CijVj, with the current flowing in wire i and 
the voltage Vj applied to the wire j. The current conservation, = 0, 
and the absence of response to the equal change in voltages result in the 
Kirchhof’s rules, — 0- suggests that we can choose more 

convenient linear combinations of Jj, Vj reducing the number of independent 
components in Cij. In the d.c. limit we have from the Kubo formula Cij = 
— Yij), with Yij = IS'ijp ; one can also write Yij = Tr(piPj). [19] 

The appropriate representation for the reduced conductance matrix may 
be constructed by using generators of 17(4) Cartan subalgebra, which are 
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Figure 1: A four-lead junction of 
quantum wires, corresponding to 
the Hamiltonian Q. 



Figure 2: Point contact of two wires 
is schematically shown, illustrating 
discrete symmetries in our choice of 
S'-matrix, Eq. (10). 


three traceless diagonal matrices and one unit matrix. We dehne 

= diag(l, -1,0, 0), /i 2 = diag(0, 0,1, -1), 

= 1/a/2 diag(l, 1, -1, -1), /i 4 = 1/a/2 diag(l, 1,1,1). 


with the property = 25^*,, j = 1,..., 4. The densities are expressed 

as pj = l/y/^J^kRjkf^k, where the 4x4 matrix R is given by 

( t . » lb 

_L 0 11 

R= w _'i ? (3) 

^ 72 2 2 

V 0 -A V 

and has the properties = R^, detH = 1. The outgoing amplitudes are 
expressed in a similar form with pj replaced by pj = S~^ ■ pj ■ S. It also means 
[20] that we work now with the combinations of currents and voltages of the 
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form I"* = or 


'new 

k 


ineu, ^ ^ yy., ^ _ yy^ ^ 

ir^ = ih - h)/2 , Vr^ = {Vs - V,)/2 , 

/3-- = (/, + J2-J3-J4)/y8, 


new 


new 


new 


new 


new 
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= {V + V2-Vs-V,)/V8, 
y j./Ts, = y Vj/Vs, 


new 


new 


new 


We could use more physical combinations Vi — V 2 etc. without additional 
factors 1/2, l/\/8, see [19], however, it is irrelevant for our purposes below. 

The reduced conductance matrix in such basis is determined by G = 
R C and has a structure 



(5) 


In the presence of interactions, the structure of the last expression is un¬ 
changed, but the elements vary. The main effect in the d.c. limit can be 
described by the renormalization of the S'-matrix, [19] which translates to 
the renormalized quantity. 



( 6 ) 


where the superscript R shows that we work in the basis pj instead of pj, 
and the superscript r denotes that the quantity is fully renormalized by 
interactions. From now on we assume that all quantities are renormalized 
and drop this latter superscript. 

4 Renormalization group equations 

The renormalization of the conductances by the interaction is determined 
by hrst calculating the correction terms in each order of perturbation theory 
in Qj. We are in particular interested in the scale-dependent contributions 
proportional to A = ln(L//), where L and I are two lengths, characterizing 
the interaction region in the wires, Eq. ([^. The above expression for A corre¬ 
sponds to vanishing temperature, T <^vf/L, and for higher temperatures we 


6 


should replace A = In(/ ^ j ^,T/vp\). In lowest order of perturbation 

the scale dependent contribution to the conductances is given by m 

Gjk ^ + ?y^Tr 

l^m 




Qml^ 1 


(7) 


where Gjfc|g=o = ^jk — Yjk, the Wjk = [Pj,Pk] are a set of sixteen 4x4 
matrices (products of hh’s are matrix products), pmi = gm^mi is the matrix 
of interaction constants and the trace operation Tr is dehned with respect to 
the 4x4 matrix space of fh’s. 

We multiply Gij by from the left and by R from the right to get the 
components of in the form 


yR _ yR\ 
jk jk\g=o 



( 8 ) 


Differentiating these results with respect to A (and then putting A = 0) we 
hnd the RG equations in the hrst order in the interaction 


jLyR 

dA 



(9) 


The number of non-zero matrices Wjl = {R^ ■ W ■ R}jfe is reduced to nine 
in most general case, since = 0 ; we have g^i = {R^ ■ g • R}m«- 

The matrices Wjl are best evaluated with the aid of computer algebra. 


5 RG flow ambiguity phenomenon 


We continue our analysis by appropriate choices of interaction gij and S- 
matrix parametrization. Let interaction constants hrst be equal in both 
wires g = dia.g{g, g, g, g), and we parametrize the S-matrix by three angles 
in the following way 


/ ri 

h 

f2 

/2 \ 

h 

n 

/2 

/2 



r2 

^2 

Ui 


G 

r2 ) 


( 10 ) 
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where 


n = +COS/?), ti = +COS/3), 

rs = + cos/?), ^2 = + cos/?), 

fi = f 2 = ^sin/?. 


( 11 ) 


This parametrization (10) describes the transmission and reflection in each 


wire (tti, 02 ) and the hopping between the wires (/?). The above Eq. ([^ 
is invariant npon “rephasing”, i.e. the mnltiplication of S from both sides 
by nnitary matrices of the form diag(e*'’'be*'’'®,e*'’''^). Withont loss of 
generality we may assnme /? G [0,7r/2], a \^2 ^ [—7r,7r]. 

The compntation of the rednced condnctances matrix yields 


_ 


^cos/?cosai 0 0 

0 cos /? cos 02 0 

0 0 cos2 /? 


oi 0 O' 

0 02 0 I (12) 

0 0 6 , 


The matrix form of the RG eqnation ([^ is now written as a set of conpled 
RG eqnations for components of in terms of the initial variables 


dA 


= g sin^ /? (cos ai(2 cos /? + cos 02 ) — 3) 
+1 — cos^ ai cos^ ~ \ sin^ /? sin ai sin 02 ) , 


dYj dYif 
dA 


dA 


(13) 


Q:i-f^a2 


33 


dA 


= sin^ 13 cos /? (l(cos ai + cos 02 ) + cos /?) . 


Onr natnral desire is to express these eqnations entirely in terms of the 
condnctances, as it was snccessfnlly done in onr previons stndies for two and 
three wires connected by the jnnction [T^l 1221 US]- 

We have three independent components of the rednced condnctances ma¬ 


trix, denoted as oi, 02 , b in Eq. (12). The attempt to write the right-hand 


side of the RG eqnations (13) in terms of oi, 02 , 6 faces the ambignity prob¬ 
lem in the term oc sinaisinQ! 2 - The sign of this term depends on angles 
range: if ai and 0:2 both belong to the range either (0, tt) or (—vr, 0) then 
sin ai sin 0:2 is positive, bnt if ai and 0:2 belong to different segments (0, tt) 
and 02 to (—vr, 0) then the discnssed term is negative. Notice, that the valnes 










of conductances Oi and 02 are not affected by the change of sign ai —)■ —ai 
and 0^2 —^ — 0 . 2-1 respectively. This change of sign corresponds to complex 
conjugation of some elements of S'-matrix (10), namely, r* —)■ r* and U t*, 
which cannot be compensated by “rephasing” operations. 

One may ask what is the internal discrete symmetry, manifesting itself 
at the level of RG equations? To answer this, we consider two decoupled 
(/? = 0) Luttinger liquids with impurities. Standard calculations [I2] show 
that (at least in the lowest Born approximation) the phase aj is equal to 
Ubs/vp, with Ubs the backscattering amplitude off the impurity in jth wire. 
The sign of Ubs is unimportant, as only the square {Ubsl"^ defines the conduc¬ 
tance [7]. When allowing hopping between the wires, we start to feel the 
difference between two cases: i) the scattering potentials in both wires are of 
the same sign, i.e. both bumps or both pits, or ii) scattering potentials in the 
wires are of different sign, i.e. a pit in one wire and a bump in another. An¬ 


other explanation of the sign question in Eq. (13) concerns the particle-hole 
symmetry of the Hamiltonian ([^. At the level of decoupled wires the sign 
of Ubs (and aj) is changed when passing to hole description, —)■ etc. 
One may then regard the above sign question as arising from the possibility 
to perform the particle-hole transformation in one wire. 

We see that the RG equations cannot be dehned in terms of conduc¬ 
tances only. Generally, we have two different RG flows for conductances, and 
the choice between them should be done on the base the initial phases of 
S'-matrix, ai, 0 : 2 • Further analysis of RG equations shows that the ambigu¬ 
ity doesn’t influence the position of hxed points (FPs). We have four FPs 
parametrized by oi = ±1, 02 = ±1, b = 1, which reads as ai ^2 = 0, vr, /3 = 0 
in terms of angles. These FPs correspond to simple cases of two separated 
wires with absolute transmission or reflection in each of them. The hfth FP 
is defined by ai = 02 = & = 0 and discussed below. 

To clarify the character of different FPs we generalize our consideration 
and allow for different interaction constants in two wires : 


g = diag(^i,^i,^2,^2). (14) 

The representation of the RG equations purely in terms of conductances is 
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cumbersome and we rewrite them in terms of the angles 


4cos/3 


(£/i(l + 3 cos^ 13) sin tti — g 2 sin^ [3 sin 02 ) 


da I 
~dA 

da 2 

dA dA 

^ ^ sin /? {gi cos ai + g 2 cos 02 + ( 5-1 + ^' 2 ) cos 13) 


dai 


( 16 ) 


310-52 


The same ambiguity of RG equations in terms of conductance is seen again 
in (15). One can change, e.g. 02 —)■ — 02 , without changing the conductance 


but altering the RG flows. 

As before, we observe four universal FPs, a \^2 = 0, tt, /3 = 0 (i.e. 01^2 = 
±1,6 = 1). Only one of these four is a stable FP and it is dehned by the sign 
of the interaction in individual wires. According to usual expectations mm 


we have for the stable FP: a,- = 


= sign gj. 


In addition we hnd a hfth non- 


universal FP, which is never stable, and whose position attains a compact 
form in terms of conductance: 


i/2 yi , 2 61 

Oi = —02 =-, 0 = a, . (_1d) 

9i + 92 

This FP is in physical domain at |ai| < 1, which happens for gig 2 > 0. 

We illustrate our hndings in Fig. where we show the body of conduc¬ 
tances, i.e. the set of allowed conductances values (oi, 02 , b) and possible RG 
trajectories, starting from the same parametric point ( 01 , 02 , 6 ) for different 
values of 91 , 92 - We see that the stable FP depends on the quarter in the plane 
of { 91 , 92 ), similar to the situation with Y-junction |2l]. We also demonstrate 
the existence of two possible RG flows, leading to the same stable FP (darker 
one is for plus sign in ambiguity term, lighter - for minus). One can verify 
that the scaling exponents near the FPs are not affected by the discussed 
sign ambiguity. Two possible RG flows result only in a different prefactors in 
the scaling dependence of the conductances. For instance, for the repulsive 
interaction in both wires, 91,92 > 0, if the RG flow starts in the vicinity of 
FP |q;i|, |a 2 |,/d -C 1 then we have 

1 - ai ,2 ~ ^ ^ 2 ^^g- 2 gi, 2 A^ 1 - 6 ~ (17) 

Starting far away from such vicinity, two possible RG trajectories will end 
at the different points aj,/3 in but the scaling law is the same. 
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Figure 3: Fixed points and RG trajectories in conductances space. Two 
possible trajectories lead to the attractive FP, determined by the signs of 
interaction in wires. It is seen that the RG flows of the conductances are 
non-monotonous functions of the scaling variable. The hfth non-universal 
FP appears when gig 2 > 0, its position dehned by Eq. (16). 


6 Discussion 

In this paper we study four-lead junction of spinless Luttinger liquid wires 
by fermionic representation in scattering state formalism. The interaction 
in wires leads to the renormalization of the conductances of the system, ex¬ 
pressed via the absolute values of the F-matrix elements. The RG equation 
for the conductances can usually be formulated entirely in terms of con- 
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ductances, which was explicitly checked for general two-lead and three-lead 
junctions. In this paper we demonstrate that in case of four-lead junctions 
the RG equations possess the sign ambiguity unresolved in terms of con¬ 
ductances. From a mathematical viewpoint, this ambiguity is an intrinsic 
property of U (4) group, and is most easily seen in the sign choice, which 
dehnes the left- and right-isoclinic subgroups of S'0(4) C U{4). From a 
physical viewpoint, the ambiguity may be traced back to the particle-hole 
symmetry of our Hamiltonian. This results in two possible non-monotonous 
dependences of the renormalized conductances as functions of the scaling 
variable. 

We note that should we use the bosonization approach in our analysis, we 
would not observe the ambiguity in question. This is because the bosoniza¬ 
tion starts with the FPs of RG equations and analyzes the scaling dimensions 
of perturbations around it. We show above that the scaling dimensions (ex¬ 
ponents) do not depend on the choice of RG flow, and it is only prefactors 
before the scaling exponentials which are determined by the particular RG 
trajectory, leading to the close vicinity of FP from the distant points in 
parameter space. 

We believe that the discussed ambiguity may be an important issue in fur¬ 
ther theoretical investigation and experimental manipulation of X-junctions 
between quantum wires. 
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